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Abstract 

Three-dimensional bond or site percolation theory on a lattice can be 
interpreted as a gauge theory in which the Wilson loops are viewed as 
counters of topological linking with random clusters. Beyond the perco- 
lation threshold large Wilson loops decay with an area law and show the 
universal shape effects due to flux tube quantum fluctuations like in ordi- 
nary confining gauge theories. Wilson loop correlators define a non-trivial 
spectrum of physical states of increasing mass and spin, like the glueballs 
of ordinary gauge theory. The crumbling of the percolating cluster when 
the length of one periodic direction decreases below a critical threshold 
accounts for the finite temperature deconfinement, which belongs to 2-D 
percolation universality class. 

1 Introduction 

Percolation is a purely geometrical phenomenon which in many respects resem- 
bles a continuous thermal phase transition. The theoretical description of the 
percolation processes is conventionally given in terms of the cluster sizes pQ, 
and most of the universal scalings deal with size distribution of clusters. The 
point of view which is taken in this paper is different. We focus on topological 
entanglement of random clusters and use it to describe how percolation theory 
in three dimensions can be viewed as a full-fledged gauge theory. 

Though the gauge group of the theory in question is trivial (it is the <7 — * 1 
limit of the symmetric group S q ), the occurrence of a confining phase yields some 
new hints on the mechanisms of quark confinement of more general theories. 
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Transcribing percolation in terms of gauge theory has also some important 
consequences for percolation itself. In three-dimensional systems there are al- 
most no exact results whatsoever, but in gauge theories we have a number of 
well- verified conjectures that can be translated into the language of percolation. 
In this way we shall, for instance, relate certain linking properties of the closed 
paths within a percolating cluster to the universal quantum fluctuations of the 
chromoelectric flux tube joining a quark pair in the confining phase of whatever 
gauge theory. 

An unsuspected property of random percolation which emerges from this 
new viewpoint is that its spectrum should be composed by a (possibly infinite) 
tower of physical states of increasing mass and spin — the glue-balls of the 
corresponding gauge theory. Their mass ratios near the percolation point define 
a totally new set of critical amplitude ratios belonging to the universality class 
of 3-D random percolation. 

Another piece of useful information comes from considering percolation in 
a slab which is infinite in two dimensions, but of finite length I and periodic 
in the remaining direction. The associated gauge model describes a system 
at finite temperature T = l/£. This transition is accurately described by the 
universality class of two-dimensional random percolation, but the corresponding 
deconfining temperature T c may be used to define a new critical amplitude of 
the three-dimensional system. 

We test the universality of this new set of critical amplitudes by performing 
large scale numerical experiments in three different kinds of lattices: the simple 
cubic (SC) with bond or site percolation and the body centered cubic lattice 
(BCC) with bond percolation. The numerical implementation of these systems 
is straightforward in comparison with simulations of ordinary gauge models: no 
Markovian process is needed and there are neither thermalization problems nor 
critical slowing down. Preliminary results have been presented in [2]. 

The organisation of the paper is as follows. In the next Section we define 
a new class of observables of the percolation theory to be identified with the 
Wilson loops of the corresponding gauge theory. In Section [31 we describe a 
method to study the linking properties of the random clusters which is at the 
heart of our analysis. Section 0] is devoted to the comparison of confinement 
mechanisms and in Section [5] we extract the string tension and show the rele- 
vance of the universal terms generated by quantum fluctuations of an underlying 
effective string. In Section El we study the plaquette correlators in order to find 
the low-lying states of the spectrum. Section is devoted to the transcription 
of the percolation on a slab into a gauge system at finite temperature and Sec- 
tion [S] discusses the universality class of the finite temperature deconfinement 
transition. Finally in Section [5] we draw some concluding remarks. 

2 Observables 

The most basic observables of any gauge theory are the Wilson loops. These 
are operators which assign to each pair (C, 7) formed by an arbitrary gauge 
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configuration G and any closed path 7 of the space a suitably defined complex 
number W 7 (G) . Their importance stems from the fact that they serve as order 
parameters for confinement: in the confining phase the vacuum expectation 
value of large Wilson loops exhibits area law. 

One can define similar observables in the framework of random percolation. 
In this context the configurations are generated simply by occupying each site 
or bond on an initially empty 3-D lattice A with independent probability p. 
Two sites are considered neighbours if they share one bond. The resulting 
configuration is a graph G drawn on the lattice, composed by the occupied 
bonds (bond percolation problem) or by the bonds joining occupied neighbour 
sites (site percolation problem). The connected components of G are the clusters 
of the configuration. We choose as 7's the closed paths of the dual lattice A. The 
value W 7 (G) measures the topological entanglement between 7 and G. More 
precisely we apply the following rule 

1. W 7 (G) = 1 if no cluster of G is linked to 7; 

2. W 7 (G) = otherwise. 

This definition did not come out of the blue. Starting from the Fortuin- 
Kasteleyn random cluster representation of the 3-D Ising model [3] combined 
with the Kramers- Wannier duality @] , it is easy to express the Wilson loop W 7 
belonging to a Z2 gauge system in terms of the winding numbers modulo 2 of 
the Fortuin-Kasteleyn clusters across the 7 loop [5]- Generalising this result 
to g-state Potts model one is led to the above rule in the case of non-integer 
q. A recently developed algorithm could even allow to evaluate explicitly these 
quantities for real q > UJ. We extend the above recipe to any percolating 
system, owing to the fact that bond percolation can be viewed as the q — > 1 
limit of the g-state Potts model. 

The removal of an occupied bond b from a graph G can lead to two different 
issues. If the number c(G) of clusters is kept invariant, then b is necessarily a 
step of a closed path (or loop) of G, whereas if c(G) increases by one, b cannot 
lie in a loop and is called a bridge (see Fig. P). Clearly only the former bonds 
contribute to W 7 (G). If two graphs G and G' have the same loops and differ 
only in the bridges, then they yield the same value of W 7 for any 7. In more 
technical terms we may write the double implication 

GUG'-GnG' is a tree <S> W 7 (G) = W 7 (G') , V7 c A , (2.1) 

hence the transformation G — * G' has some resemblance to a gauge transforma- 
tion. 

Note that the connected correlator among occupied bonds is exactly zero 
by construction, but this is not an invariant quantity under the above G — ► G' 
transformation. Cutting all bridges of G yields the maximal invariant subset 
Bq of G, made by bonds belonging to some loop. Of course we have W 7 (G) = 
W 7 (Bq) for any 7. The connected correlator among bonds belonging to Bq is 
by no means trivial and is directly related to the connected correlator of the 
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W=l 




Figure 1: The closed thick lines 7 and 7' represent Wilson loops. The dashed 
lines are bridges of the cluster. The other solid lines are closed paths of the 
cluster. 7 is linked to the cluster, while 7' is unlinked. 



plaquette, i.e. the Wilson loop Wu associated to the smallest loop □ C A, 
because Wrj = if and only if the bond dual to □ belongs to Bq. In Section 
we shall use such a correlator to extract the low-lying mass spectrum of the 
theory. 

3 Cutting all bridges 

In our approach the only bonds which play a role in the evaluation of the 
observables defining the gauge theory are those belonging to loops. Thus, once 
a new configuration G is generated, we first get rid of all bridges. One way of 
achieving this goal is the following. 

1. Eliminate all the dangling ends (see Fig. [5]} and b). At this stage the 
remaining graph is formed by loops and lines of bridges connecting them. 

2. Build a reduced graph in which the only vertices are the lattice sites with 
more than two incident occupied bonds. The edges of the graph are formed 
by lines of occupied bonds (Fig. EJd). 

3. Erase one edge at a time and apply each time a cluster reconstruction 
algorithm (for instance Hoshen-Kopelman [7]) in order to pick off the 
remaining bridges (Fig. [SJi). 

In order to check whether a planar loop 7 C A is linked to a configuration 
G, we first project out all bridges as discussed. Then we switch off the layer 
of occupied bonds which pierce the planar surface £ encircled by 7 and rebuild 
the cluster structure of the cut graph. For a non-trivial linking there must be 
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a b c 

Figure 2: Cutting all bridges of a configuration 




Figure 3: A 2-D sketch of the method used to evaluate the linking of a configu- 
ration with a loop 7. The two lines of solid dots represent the sites on either side 
of the layer of bonds piercing the surface E (dotted line). The vertical broken 
lines are the switched off bonds. The big open circles represent the clusters of 
the cut graph to which the sites of the layer belong. The cluster a reaches both 
margins of the layer. 



at least one cluster which reaches the layer on either side (see the graph on the 
left of Fig. |3J) ■ In such a case we build an auxiliary graph in which the vertices 
represent the clusters of sites on either side of the layer; a cluster connecting 
sites of opposite sides is represented by two vertices, one for each side (see the 
graph on the right of Fig. |3J) • We draw an edge between two clusters in the 
opposite sides of the layer if they are connected by switched off bonds. The 
configuration is truly linked to 7 if there is a path (at least) joining two vertices 
lying on opposite margins of the layer, but belonging to the same cluster of the 
cut graph (like the vertices a and a' of the figure). 

The whole procedure of cutting all bridges of a configuration and then eval- 
uating its linking property with a set of Wilson loops is time demanding, thus 
a good implementation is mandatory. However the variance of the measured 
quantities turns out to be small, hence one can reach more precise results than 
the corresponding estimates in ordinary gauge theories. 
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4 Confinement mechanisms 



According to the recipe given in the previous Section, the vacuum expectation 
value of the Wilson operator W 1 is defined as 

(W 7 ) = lim Y W,(G t )/N = number of config. unlinked to 7 ^ 
jv-»oo total number 01 configurations 

What is the functional form of this quantity for large loops? It depends on the 
value of the occupation probability p. If p is below the percolation threshold 
p c , then there are only finite-size clusters. If the loop 7 has much a larger size 
than the clusters, then the configurations where W 1 (Gi) = ncccsarily have 
some cluster located near the loop. The number of these clusters grows linearly 
with the perimeter I7I of 7 and produces the exponential decay (W 7 ) oc e~"' 7 '. 
We say that the theory is deconfincd. On the contrary if p > p c the theory is 
confined. Confinement is expected to show up in an area law for the expectation 
value of large Wilson loops. Indeed in such a case there is an infinite cluster and 
the number of closed paths linked to 7 grows with the area A of the minimal 
surface £ encircled by 7. These paths will pierce £ at points; let N = aA be 
their mean number. Assuming these points to be randomly distributed on the 
surface, the probability of finding k such points inside £ is binomial, 

P N (k)=( I l^a k (l-a) N - k . (4.2) 

Note that only the k = term contributes to the numerator of (|4.1(l . so the 
expectation value of the Wilson loop becomes 

(W 7 ) = (1 - a) N = e- aA , a = -a log(l - a) . (4.3) 

One thus apparently obtains an area law decay with string tension a for any 
Wilson loop, including those of small size. There is however a flaw in the 
argument; even if the configurations are obtained by populating each bond (or 
site) of the lattice independently with a probability p, when all the bridges are 
erased it is no longer true that the remaining bonds are randomly distributed, 
as anticipated previously. In fact, since the interaction among the intersection 
points of £ is rather weak, we expect an area law only for large enough Wilson 
loops. 

In gauge theories two different confining mechanisms were proposed. One 
is based on the condensation of center vortices [Hj- These objects are string- 
like structures which are created by gauge transformations with a non-trivial 
homotopy associated to the center of the gauge group C(G). In a 3-D lattice the 
center vortices are represented by a skein of loops forming an infinite network in 
the confining phase |5J . Each center vortex linked with a Wilson loop contributes 
to it with a factor of £ € C(G), thus each configuration contributes with a factor 
C™, where n is the number of linked loops. Although at the end the overall effect 
is again a decay of (Wj) with an area law, such a mechanism is slightly different 



6 



from the one we have described in pure percolation, where a single linked loop 
suffices to gives a weight zero to the configuration. Another difference is that 
center vortices carry some conserved charge. For instance, if C(G) = Zjy the 
vortex flux is conserved modulo N, while paths of random percolation do not 
carry any conserved charge and can intersect freely. 

The other confining mechanism is based on the old conjecture JD] that the 
vacuum behaves like a dual superconductor. The key element of this picture is 
the monopole condensate which squeezes the gauge field generated by a pair of 
sources (quarks) into a thin flux tube (the dual version of the Abrikosov vortex). 
This causes the Wilson loop to decay with an area law. One is led to conjecture 
that such a thin flux tube should vibrate as a free string ^Tj. As a consequence, 
the expectation value of a rectangular Wilson loop of size R x T is expected to 
have the following asymptotic functional form in the continuum limit |f 2| 

(W(R, r)) = C e -K* + T)-^T f (4.4) 

where C, p and a arc functions of the coupling constant and 77 is the Dedekind 
function 

00 

V(r)=q 1/2i l[(l-q n ) , q = e*™ . (4.5) 

n=l 

The factor under the square root accounts for the universal quantum contribu- 
tions of the supposedly string-like flux tube describing the interaction between 
far-away sources. 

Random percolation, lacking any non-trivial conserved charge, can hardly 
account for effects which play the role of magnetic monopole condensation. 
Notwithstanding this difficulty, we get indirect evidence of the formation of 
a vibrating string-like flux tube by measuring the universal shape effects it 
produces, as discussed in the next Section. 

5 String tension 




0.0001 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

5 10 15 20 25 30 35 

Figure 2] Square Wilson loop as a 
function of R for bond percolation in a 
cubic lattice of size L z = 64 3 at 
p = 0.26 . 



0.2 0.4 0.6 

Figure^] The universal shape effects of 
Eq. 15.31 are compared with numerical 
data of Eq. 15.21 for three different 
values of p. 
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We estimated the string tension a by fitting the mean values of the Wilson 
loops associated to squares of side R to Eq. 14.4(1 . that in such a case becomes 



W(R) = C R* exp(-2p R — a R 2 ) 



(5.1) 



Typically, in a lattice of size L 3 we considered all the squares with R < L/2. The 
fits for not too small R are very good (see Fig. 0}, nevertheless the parameters 
slightly depend on the value R m in of the smallest square included in the fit. 
Since these formulae are expected to be valid only asymptotically for large values 
of R, we progressively eliminated the data of lower R until stable parameters 
were obtained. 

In order to check for presence of universal shape effects ascribed to the quan- 
tum fluctuations of the effective string, we considered, as in |13| . the quantity 



K(n,R) = e 



.W{R + n,R-n) 



(5.2) 



which asymptotically (i.e. large R and R 
ratio t = ^, namely, 

K(n,R)^f(t) = * 



W{R) 

- n) should be only a function of the 



(5.3) 



Note that it does not contain any adjustable parameters. This function is 
plotted in Figure and compared with the numerical data for three different 
values of p. The presence of the expected universal shape effects seems uncon- 
troversial. 

The string tension a is a physical quantity with the dimensions of an inverse 
square length, hence it is expected to exhibit the following power law behaviour 
sufficiently close to the percolation threshold 



S(j>-Pc) 



2i> 



(5.4) 



where v is the critical exponent of the correlation length in 3-D percolation. We 
used the value v = 0.8765(16)(2) of Ref. [H] 1 . 



Lattice 


Pc 


S 


X 2 /d.o.f. 


SC site 


0.3116081(7)(2) [14 J2 


3.370(8) 


1.15 


SC bond 


0.2488126(5) 15 


8.90(3) 


0.30 


BCC bond 


0.1802875(10) 15_ 


22.07(2) 


0.98 



Table 1: The amplitude of the string tension for three different lattices. Errors 
in parenthesis affect the last digits. 

1 The first number between parenthesis is the statistical error, the second comes from the 
uncertainty in the scaling correction exponent u). 




Figure 6: The string tension a as a function of p — p c for three different lattices: 
BCC bond (top line) SC bond (middle line) and SC site. The three parallel lines 
are one-parameter fits to Eq. 115.4(1 . The corresponding amplitudes are reported 
in Tabled 



The percolation threshold p c depends on the lattice and on the kind of per- 
colation (site or bond) which is studied. We checked Eq. (|5.4I) in the percolating 
region of three different lattices: simple cubic (SC) with site or bond percola- 
tion and body-centered cubic (BCC) with bond percolation. Precise estimates 
of p c are known and are reported in Tab. d I n the SC cases we used a lattice 
of size L 3 — 64 3 , while in the BCC case we had L 3 — 55 3 . In all cases the 
one-parameter fit to Eq. (|5.4|) is very good and the scaling window seems rather 
wide (see Fig.[BJ). The resulting amplitudes S for the three lattices are reported 
in Tab. □ 

Noteworthy, in order to extract the string tension a we have thus far assumed 
that the square Wilson loops obey the asymptotic form 1)5. lj) . where the factor 
i?* accounts for the contribution of the string fluctuations. If one neglected this 
factor and only took the area term into account, the x 2 test of the critical power 
law (|5.4|l would grow worse by one or two orders of magnitude, depending on 
the kind of lattice. We consider this fact another strong evidence of a vibrating, 
confining string. 
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6 Spectrum 

As for standard gauge theories, we expect that the confining phase of percolation 
possesses a rich spectrum of physical states with increasing mass and spin, that 
we still call glueballs. 




1e-06 1 1 1 1 1 1 1 

2 4 6 8 10 12 

R/a 

Figure 7: Comparison between the standard cluster correlator defined in Eq. 
1)6.5(1 and the connected plaquette-plaquette correlator in a 32 x 32 x 64 cubic 
lattice at p = 0.260. The straight line is an exponential fit to the cluster 
correlator data (x). These data have been displaced downwards by two orders 
of magnitude for clarity. The other line is a two-exponential fit to the plaquettc 
data (+). The semi- logarithmic plot makes it evident that in the latter case a 
single exponential does not suffice. 



The basic method that goes into the computation of such mass spectrum is 
very simple. One first constructs a linear combination of Wilson loops on a fixed 
time slice of the three-dimensional lattice carrying the quantum numbers of the 
state one wishes to investigate. One then builds zero momentum operators by 
summing such a linear combination over the entire spatial lattice. The simplest 
example of a zero momentum operator coupling to the spin + states is given 

by 

$ 0+ W=^[^ni(z,2/,*) + *W*,J/,i)] , (6.1) 

x.y 

where W\jj(x) denotes an elementary plaquette variable with base at x = 
{x,y,t) and orthogonal to the j coordinate axis. According to Section |2 in 
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random percolation the plaquette variable is replaced by the dual link vari- 
able ij(x), defined as equal to 1 if the corresponding bond belongs to the 
subset Bq of loops and null otherwise. The low-lying mass spectrum can 
be extracted by studying the exponential decay of the connected correlator 
C(t) = ($ 0+ (t)<f> 0+ (0))-(<I> 0+ ) 2 , which is expected to have the following asymp- 
totic expansion 

C{t)=Y,c n e- m - t , (6.2) 

n 

where c n denote positive constants and m n are the glueball masses. An example 
of such a correlator in a SC bond lattice is reported in Fig. [7] where it is evident 
that at least two different scalar states contribute to C(t). The estimates of the 
lowest mass in the range 0.258 < p < 0.270 fit well with the scaling form (see 
Fig.© 

m = M (p- Pc r , (6.3) 

with 

M = 12.45 ± 0.07 . (6.4) 



0.45 



0.35 



09 
09 

a 



0.25 



0.15 



0.256 



0.26 



0.264 



0.268 



0.272 



Figure 8: The mass of the lowest state as a function of p in a simple cubic bond 
percolating lattice. The dashed curve is a one-parameter fit to Eq. (|6.3|l . 



Also the first excited state seems to follow the same power law as expected, 
though the errors are rather big. Its mass Mq is about twice M . 

Such a behaviour is very different from the one observed in the standard two- 
point correlation function of the percolation problem, defined as the probability 
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G(x,x') that the sites x and x 1 are in the same cluster 2 . The corresponding 
connected, zero-momentum, projection 



C(i) = £ [G(0,jg)-G(0, oo)] 



(6.5) 



exhibits a single exponential behaviour with a mass term which coincides, within 
the numerical accuracy, with that of the lowest energy state coupled to the pla- 
quette operator (see Fig.[7J). This indicates that the standard cluster correlator 
couples only to the lowest glueball, while the new observables suggested by 
the gauge theory interpretation of the percolation disclose a totally unexpected 
spectrum of physical states. 

The lightest spinning glueball is the 2 + state. It can be observed in the 
exponential decay of the correlation function of the operator 



A difficulty encountered in this case is that the signal is drowned within the 
statistical noise for values of t beyond three or four lattice spacings. In spite of 
this accuracy problem, one can still verify that in a SC bond lattice the scaling 
form l|6.3[) is approximately obeyed with an amplitude 



In order to have more accurate results on the mass spectrum the basis of the op- 
erators should be enlarged to Wilson loops of different shapes, trying to enhance 
their overlap with the glueball states. 

7 Deconfinement at finite temperature 

In quantum field theory the concept of temperature is introduced by simply 
compactifying the Euclidean time direction and identifying the inverse temper- 
ature with the temporal extension of the space-time manifold. Lattice field 
theories at a temperature T are formulated in a slab which is infinite in the 
spatial dimensions, but of finite length I = 1/T (in lattice spacing units) and 
periodic in the remaining temporal direction. 

In any confining gauge theory there is a critical temperature T c above which 
the system is deconfined in the sense that for T > T c the string tension a van- 
ishes. In this Section we demonstrate that the same phenomenon also occurs 
in random percolation. In the latter case the deconfinement mechanism is par- 
ticularly transparent, showing its purely geometric origin: Non- vanishing string 
tension requires an infinite, percolating, cluster. Shrinking the width of the slab 
reduces the number of possible percolating paths along the spatial directions 
until the infinite cluster crumbles away, yielding a vanishing <j. 

2 This quantity is directly related to the correlator in q — > 1 limit of the g-state Potts model, 
see for instance Ref. 1 1 01 . page 156. 




(6.6) 



M 2 + = 80 ± 10 . 



(6.7) 



12 



To put it in different terms, note that as the temperature varies from zero to 
infinity a three-dimensional system is gradually dominated by two-dimensional 
behaviour; in particular the percolation threshold is a decreasing function of the 
space dimensions. For instance in the SC bond lattice at T = the percolation 
threshold is at p c (D = 3) = 0.2488... (see Table [TJ. At T = oo the system 
reduces to a square lattice, where p c {D = 2) = |, hence heating a system 
which at zero temperature lies in the percolating phase with p < | inevitably 
undergoes a deconfincmcnt transition at a finite temperature T c . 



1 



0.5 



0.335 0.345 0.355 0.365 

P 

Figure 9: The wrapping probability Rl in a simple cubic site percolating lattice 
for I = 7 and L = 60, 70, 100. The vertical arrow denotes the estimated value 
of p£ at L = oo and the horizontal line indicates the exact planar value of Roo . 



In order to estimate T c in various site or bond percolation lattices we con- 
sidered a slab of size L x x L y x i with L x = L y = L 3> I and periodic boundary 
conditions in all directions (see Fig. W^ . We calculated the probability Rl(p) 
for a cluster to wrap around one of the large dimensions. Wrapping probabilities 
can be defined in different ways: wrapping around the x direction, around the 
y direction, around cither direction, around both directions, etc. To be definite, 
we consider wrapping around the x direction. For large L it coincides with the 
probability that the system percolates along x. 

We evaluated Rl(p) using a very efficient algorithm described by Newman 
and Ziff ^7] . For a bond percolation problem it consists of repeatedly adding a 
random bond to an initially empty lattice, identifying the clusters joined by the 
bond and merging them if they are different. At each step one checks whether 
the touched cluster wraps around x using a clever method described in Ref. j!8| . 



13 



The process stops as soon as a wrapped cluster is detected. In this way one can 
evaluate the probability Qi(n) that a random configuration with n occupied 
bonds is wrapping around L x for any n < N, where N is the total number 
of bonds of the lattice. This method may be adapted to site percolation in a 
straightforward way. Then one simply finds the required quantity Rl (p) for any 
value of p by convolution with the binomial distribution 

Figurc[5]shows some examples of Rl(p) in the case of SC site percolating lattice. 
In the L — > oo limit this wrapping probability becomes a step function. We have 

Ro (p) = {°, l° TP<Pe , (7.2) 
\l for p>p e ' v ' 

where the threshold value p = pi depends on the type of lattice and on its 
width £. When L = oo the slab system is equivalent to a 2-D torus, where the 
wrapping probability at criticality has been calculated exactly by Pinson |19| 
and is Roo{pi) = 0.521058290... . We can use his result to measure the value of 
Pi in a slab of width £ by finding the value of p for which 

Rl(p) = Rooipt) ■ (7.3) 

This method was first applied in the case of the Ising model in Ref. [2"U| 
(with the appropriate value of Roo, of course). These estimates in the case of 2- 
D percolation turn out to scale particularly well with the system size: Newman 
and Ziff argued that in planar lattices the leading order finite size correction 
goes like 

p = Pe + cL- 2 - 1 ^ = Pe + cL- u / 4 , (7.4) 

where v-i = | is the thermal exponent of 2-D random percolation. 

We checked it in the slab geometry finding good agreement for L large 
enough, as Figure ITU1 shows in the case of a slab of width £ = 6. 

We estimated in this way the threshold pi for seven different lattices. The 
results are reported in Tab. [3 The slab widths (or equivalently the inverse 
temperatures) were chosen in such a way that the values of pi lie in the scaling 
region of the string tension tr, as determined in Section [S] In this manner we 
were able to evaluate also the ratio T c j y/a using the amplitudes of Tab. It 
turns out that these ratios for different lattices and widths are clearly compatible 
with a common value, as required by universality (see last column of Tab.|2J. 

In (d + l)-dimensional field theory at finite temperature there is a charac- 
teristic interplay between d + 1 and d critical behaviours. This is particularly 
evident in the present instance: the pi values are extracted through a 2-D per- 
colation power law l|7.4() in order to take into account the finite size scaling tied 
to L. However the tower of pi values as a function of the width I obeys a typical 
power law of 3-D percolation: 
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Lattice 


1/T 


Pi 




SC site 


7 


0.3459514(12) 


1.494(11) 


BCC bond 


3 


0.21113018(38) 


1.497(10) 


BCC bond 


4 


0.20235168(59) 


1.506(11) 


SC bond 


5 


0.278102(5) 


1.480(12) 


SC bond 


6 


0.272380(2) 


1.492(13) 


SC bond 


7 


0.268459(1) 


1.500(13) 


SC bond 


8 


0.265615(5) 


1.504(14) 



Table 2: The critical pi for different lattices at different temperatures and the 
corresponding universal ratio T c / \fij as obtained by combining Eq.s l|5.4|) and 
H7.5J) . Errors in parenthesis affect the last digits. 



where v indicates, as in all the other formulae of this paper, the thermal expo- 
nent of 3-D percolation, p c is the critical threshold as listed in Tab. ^ an d the 
amplitude T c depends on the kind of lattice. Figure ^] shows the four pg values 
of the bond SC lattice as a function of l~ x ' v ' . A one-parameter fit to Eq. (|7.5fl 
yields 

T c = 4.45 ±0.02 . (7.6) 

It should be noted parenthetically that, as thep/s are essentially two-dimension- 
al quantities, they can be evaluated with high precision. Comparison between 
Table Q] and Table [21 shows that in some cases the level of precision of those 
overcomes that of the best estimates of p c . Moreover, with a modest additional 
computational effort, it would be possible to further improve their precision. 
Thus one could perhaps envisage to apply systematically Eq. 1|7.5[) to improve 
the estimates of p c and/or v. 





t igurc tlOl Finite -size scaling of the 
estimated L — > oo value of pi for bond 
percolation in a slice of SC lattice of 
size L x L x 6 with 32 < L < 80. Each 
point corresponds to about 10 s 
configurations. The solid line is a fit to 
Eq. (fT4jl . 



b igure mi Finite -size scaling of pe as a 
function of l -1 /", where i is the width 
of a SC lattice in bond percolation. The 
solid line is a one-parameter fit to 
Eq. 17. 5H . The statistical errors are 
much smaller than the symbol size. 
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8 Universality class of deconfinement 



The deconfincd phase of any (d+ l)-dimcnsional gauge theory at finite temper- 
ature is characterised by the vanishing of the string tension a. The interaction 
between static sources (quarks) is described in terms of Polyakov operators. 
These are straight Wilson loops wrapped around the short periodic direction 
I (see Fig. EJ. The Polyakov-Polyakov correlator (P(0)P(R)) of two parallel 
Polyakov operators only depends on their relative positions in the d-dimensional 
sub-lattice. At the deconfining point it is expected to obey a power law dictated 
by the universality class of the transition. 

The critical behaviour of gauge theories at the deconfining temperature is 
well described by the Svetitsky-Yaffe (SY) conjecture [2J which can be for- 
mulated as follows. Suppose a d + 1 dimensional gauge theory with gauge 
group G has a second-order deconfinement transition at a certain temperature 
T c ; then its universality class is the same of the order-disorder transition of a 
d-dimensional spin system with a global symmetry group coinciding with the 
center of the gauge group. In particular, the Polyakov-Polyakov correlator cor- 
responds to the spin-spin correlator of such a <i-dimensional system. Therefore, 
at the critical point, it should decay as 

, , , , -* v , const. 

(p(o)p(R)) = -^r n , (8.1) 

where r\ is the magnetic exponent of the spin model. The validity of this con- 
jecture has been well established in a large number of numerical studies. 

In the present case the above conjecture requires some generalisation, ow- 
ing to the fact that not only the center, but the whole gauge group is trivial. 
Actually the SY conjecture is somehow related to the dimensional crossover in 
a layered lattice system [221 - The universality class of such a system, as it ap- 
proaches a critical point, depends on the number of spatial directions which are 
going to infinity in the thermodynamic limit. This simple observation provides 
the basis for arguing that the critical behaviour of percolation in a slab of finite 
thickness is well described by the 2-D percolation universality class. This is also 
supported by the fact that finite size scaling of threshold probability pi is driven 
by the thermal exponent v-i-, as Eq. (|7.4|) and Figure ITUI show. 

Such a conclusion is much less obvious, and more interesting, when consid- 
ering the Polyakov-Polyakov correlator. This quantity can be evaluated using 
exactly the same method described in Section |21 for the Wilson loops. One sim- 
ply has to take into account that in the Wilson loop the surface £ (see Fig. 0J 
used to evaluate cluster wrapping is a rectangle, while in the latter case is a 
cylinder bounded by the the two Polyakov lines. Periodic boundary conditions 
put into play another difference: there are now four topologically different sur- 
faces bounded by the two Polyakov lines (see Figure EH, hence the correlator 
must be written as the sum of these four contributions. When the distance R 
between Polyakov lines is much smaller than lattice size L the main contribution 
comes from the top two surfaces of Figure ED 
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One instance is reported in Figure^] where we plot the estimated Polyakov- 
Polyakov correlator extracted from 10 4 configurations in a 200 x 200 x 6 lattice 
at the critical value of pq as determined in Tab. [21 Plotting these data versus 
R~ n with r\ — shows linear behaviour, as expected for a critical system which 
lies in the universality class of 2-D percolation. 




L 

Figure WI\ The slab geometry of the 
finite temperature setting. The two 
thick lines represent two Polyakov loops 
wrapping in the periodic temperature 
direction. 



Figure ED The dashed lines represent 
the 1-D section of the four topologically 
different surfaces bounded by the two 
Polyakov loops (solid circles) in the case 
of periodic be. 



9 Conclusions 

We have studied some consequences of a new point of view in three-dimensional 
random percolation which allows us to reinterpret it as a full-fledged gauge 
theory. 

A key difference between the conventional and the present approach is that 
instead of studying the universal scalings through the size distribution of random 
clusters we only consider their topological entanglement with suitable closed 
paths. This suggests a reformulation in the language of gauge theory. In this 
context a more detailed description of 3-D percolation universality class emerges 
as a major issue: the gauge-percolation dictionary leads to define new physical 
quantities that can be used to extract new universal amplitude ratios. 

A typical example is the ratio T^/cr between the (square of) deconfining 
temperature T c and the string tension a. We evaluated it on seven different 
lattices (see Tab. |5J finding excellent agreement with universality. Notice that 
neither T c nor a are truly foreign concepts of percolation theory: slab perco- 
lation has always been a subject of intensive study. Strictly speaking, string 
tension does not seem to have been considered previously in percolation stud- 
ies, being a typical notion of gauge theory. There is however an intimately 
related quantity, the surface tension, which can also be defined in percolation 

■ The novelty introduced by the gauge theory interpretation is a previously 
unsuspected relationship between slab percolation and surface tension. 
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0.65 ' — 
0.6 



R 

Figure 14: Polyakov-Polyakov correlator at deconfing temperature as a function 
of the scaling variable R~ v , where rj = ^ is the magnetic exponent of 2-D 
percolation. The lattice size is 200 x 200 x 6. The solid line is a one-parameter 
fit to Eq. (H3J of the data with R > 2. The resulting x 2 /d.o.f. is 0.46. 



Another class of universal amplitude ratios came to us as a surprise. It 
turns out that the Wilson loop correlators receive contribution from a tower of 
physical states of increasing mass and spin, like in ordinary gauge theories. The 
ratios of their masses define universal quantities which further characterise the 
universality class of 3-D percolation. 

Transcribing percolation into gauge theory language has also some interest- 
ing consequences in the study of quark confinement mechanisms. The confine- 
ment generated by a random percolating cluster is similar, but not identical, 
to that produced by an infinite network of center vortices. In particular the 
latter carry some conserved charge which regulates their mutual intersections, 
while the former do not carry any conserved charge and intersect freely. This 
indicates that the intersection rules of center vortices do not play an important 
role in producing confinement. 

One of the most surprising findings of the present approach is the observation 
of some universal shape effects in the vacuum expectation value of Wilson loops 
which have been always ascribed to a different picture of confinement. This 
picture says that the flux of the gauge field generated by pair of quarks is 
squeezed by the magnetic monopolc condensate into a string-like structure which 
can vibrate freely. It turns out that these very vibrations generate the universal 
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effects mentioned above. Apparently, the two different pictures of confinement 
are different descriptions of the same physical phenomenon. 

It would be interesting to extend our percolation approach to a 4-D gauge 
theory. This would require the study of plaquettc percolation and the Wilson 
loops should measure their entanglement with closed surfaces. 
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